Generalised isothermal models with strange equation of state 



X 
J3 



S. D. Maharaj and S. Thirukkanesh^ 

Astrophysics and Cosmology Research Unit, 
School of Mathematical Sciences, 
University of KwaZulu-Natal, 
Private Bag X54001, 



G^ ■ Durban 4000, 

Q ■ South Africa. 

(N ■ 

D 

Abstract 



We consider the linear equation of state for matter distributions that may be applied to 
qh| strange stars with quark matter. In our general approach the compact relativistic body al- 

' lows for anisotropic pressures in the presence of the electromagnetic field. New exact solutions 

are found to the Einstein-Maxwell system. A particular case is shown to be regular at the 
stellar centre. In the isotropic limit we regain the general relativistic isothermal universe. We 
. show that the mass corresponds to values obtained previously for quark stars when anisotropy 

■ and charge are present. 
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1 Introduction 



In an early and seminal treatment the existence of quark matter in a stellar configuration 
in hydrostatic equilibrium was suggested by Itohi. Subsequently the analysis of strange 
stars consisting of quark matter has been considered in a number of investigations. Strange 
stars are likely to form in the period of collapse of the core regions of a massive star after a 
supernova explosion which was pointed out by Cheng et al^. The core of a neutron star or 
proto-neutron star is a suitable environment for conventional barotropic matter to convert 
into strange quark matter. Regions of low temperatures and sufficiently high temperatures are 
required for a first or second order phase transition which results in deconfined quark matter. 
Another possibility suggested by Cheng and Dai- to explain the formation of a strange star 
is the accretion of sufficient mass in a rapidly spinning dense star in X-ray binaries which 
undergoes a phase transition. The behaviour of matter at ultrahigh densities for quark matter 
is not well understood: in an attempt to study the physics researchers normally restrict their 
attention to the MIT bag model (see the treatments of Chodos et al^, Farhi and Jaffe^ and 
Witten— ). The strange matter equation of state is taken to be 



where p is the energy density, p is the pressure and B is the bag constant. The vacuum 
pressure B is the bag model equilibrates the pressure and stabilises the system; the constant 
B determines the quark confinement. The studies of Bombaci— , Li et al^^, Dey et al—, 
Xu et ai— li^, Pons et al— and Usov^^ directed at particular compact astronomical objects 
suggest that these could be strange stars composed of quark matter with equation of state 



Mak and Harko*^ found an exact general relativistic model of a quark star that admits 
a conformal Killing vector. This was shown by Komathiraj and Maharaj— to be a part of a 
more general class of exact analytical models in the presence of the electromagnetic field with 
isotropic pressures. The role of anisotropy was investigated by Lobo^, Mak and Harkoi^ 
and Sharma and Maharaj-2fi for strange stars with quark matter with neutral anisotropic 
distributions. It is our intention to study the Einstein-Maxwell system with a linear equation 
of state with anisotropic pressures; this treatment would be applicable to a strange stars which 
are charged and anisotropic which is the most general case. In §2, we write the Einstein- 
Maxwell system in an equivalent form using a coordinate transformation. A new exact 
solution, in terms of simple elementary functions, is given in §3. In addition, we demonstrate 
that it is possible to find a particular model which is nonsingular at the stellar origin. The 
limit of vanishing anisotropy is studied in §4 and we regain the isothermal universes studied 
previously. In §5, we consider the physical features of the new solutions, plot the matter 
variables for particular parameter values and show that the quark star mass is consistent 
with earlier treatments. Some concluding remarks are made in §6. 
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2 Basic equations 



It is our intention to model the interior of a dense realistic star with a general matter distri- 
bution. On physical grounds we can take the gravitational field to be static and spherically 
symmetric. Consequently, we assume that the gravitational field of the stellar interior is 
represented by the line element 

ds^ = -e^-i-)dt' + e^^^'Ur^ + r\de'' + sin^ OdcP^) (2) 

in Schwarzschild coordinates (x") = {t,r,6,(j)). We consider the general case of a matter 
distribution with both anisotropy and charge. Therefore we take the energy momentum 
tensor for the interior to be an anisotropic charged imperfect fluid; this is represented by the 
form 

Tij = diag (^-p - ]^E\pr - ^E\pt + \E\pt + \e^^ , (3) 

where p is the energy density, pr is the radial pressure, pt is the tangential pressure and E is 
the electric field intensity. These physical quantities are measured relative to the comoving 
fluid velocity = e~'^5Q. The line element ([2]) and the imperfect matter distribution ([3]) 
generate the Einstein field equations; the field equations can be written in the form 



1 r 



r(l-e-2^)]' = p+^E^ (4) 
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:^(l-e-2^) + — e~2^ = Pr-'-E\ (5) 
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lie-\r^Et (7) 



where primes denote differentiation with respect to r and a is the proper charge density. 
We are utilising units where the coupling constant = 1 and the speed of light c = 1. 
The Einstein-Maxwell system of equations (H])-© describes the gravitational behaviour for 
an anisotropic charged imperfect fluid. For matter distributions with p^ = pt (isotropic 
pressures) and E = (no charge) we regain Einstein's equations for an uncharged perfect 
fluid from (gl)-©. 

An equivalent form of the field equations is obtained if we introduce new variables: the 
independent variable x and new functions y and Z. These are given by 

X = Cr\ Z{x) = e-2^W and A^y\x) = e^^^, (8) 

which was earlier used by Durgapal and Banner] i— to describe neutron stars. In terms of 
the new variables the line element ([2]) becomes 

ds^ = -A'^y^dt'^ + -r^dx^ + ^{de^ + sin2 Od^"^), (9) 
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where A and C are arbitrary constants. The transformation ([8]) simpUfies the field equations, 
and we find that the system (H])-© can be written as 
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where dots denote differentiation with respect to the variable x. 
The definition 
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(10) 

(11) 
(12) 

(13) 
(14) 



represents the mass contained within a radius r which is a useful physical quantity. The mass 
function (I14D has the form 



1 



m[x) 



\fwp{w)d'w, 



(15) 



4C3/2 7o 

in terms of the new variables introduced in ([5]). 

On physical grounds we expect that the matter distribution for realitic stellar matter 
should satisfy a barotropic equation of state pr = Pr{p)- For the investigations in this paper 
we assume the particular equation of state 



Pr 



ap — f3, 



(16) 



where a and /? are constants. This is a simple linear relationship with desirable physical 
features and contains models investigated previously. Now it is possible to rewrite (jl0p - (jl3p 
as the system 
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(17) 
(18) 
(19) 

(20) 
(21) 
(22) 



where the quantity A = pt—Pr is defined as the measure of anisotropy. The Einstein-Maxwell 
equations as expressed in (jl7p -(|22p is a system of six nonlinear equations in terms of eight 
variables {p,Pr,Pt, A, E, a, y, Z). The system (fT7|) - (j22|) is under-determined so that there are 
different ways in which we can proceed with the integration process. Here we show that it 
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is possible to specify two of the quantities and generate an ordinary differential equation in 
only one dependent variable in the integration process. This helps to produce a particular 
exact model. 



3 New solutions 



In this paper, we choose physically reasonable forms for the gravitational potential Z and 
electric field intensity E and then integrate the system (|17 p - (|22p to generate exact models. 
We make the specific choices 

1 



^2 



a + bx"' 
2k{d + 2x) 



(23) 
(24) 



a + bx'^ 

where a, b, d, n and k are real constants. The potential Z is regular at the origin and con- 
tinuous in the stellar interior for a wide range of values for the parameters a, b and n. The 
electric field intensity E \sa. bounded and decreasing function from the origin to the surface of 
the sphere. Therefore the forms chosen in (j23p - (j24p are physically acceptable. These specific 
choices for Z and E simplify the integration process. Equation ()2ip can be written as 

a bnx"-^^ 
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(25) 



where we have used (j23|) and (p4|) . This has the advantage of being a first order linear 
equation in the gravitational potential y. 

Equation (|25p can be integrated in closed form to give 
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and I? is a constant of integration. Now from (|23p. (|24p and (j26p we can generate an exact 
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model for the system pT|) - (p2]) as follows 



e^^ = a + 6x", (27) 



e 
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^2^2^ t''-^>f ^ ^^ny exp[2F(x)], (28) 

/9 _ (a - 1) + 26n3;"-i k{d + 2x) 
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= ap- P, (30) 

pt = + A, (31) 
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The equations (|27|) - (|33| ) represent an exact solution to the Einstein-Maxwell system (|17p - (|22p 
for a charged imperfect fluid with the linear equation of state pr = ap — (3. Exact solutions 
with the equation of state pr = ap — f3 have been used to model compact objects such as 
strange stars, as shown by Sharma and Maharaj^*^, and dark energy stars which are stable, as 
demonstrated by Lobo^^. The exact solution ([27|) - ([33]) may be regarded as a generalisation 
of an isothermal universe model as we indicate in §4. 

The solution (j27p - (|33p admits singularities at the stellar centre in general. The singularity 
may be avoided for particular parameter values. If we set a = 1 and n = 1 then we generate 
the line element 

ds^ = -A^D^{1 + 5x)" exp[2F{x)]dt^ + (1 + bx)dr^ + r^{d9^ + sin^ 

where F{x) = — ^(2 + bx) + ^^4"'' [bx — k{dx + x^)], with energy density 

p 6(3 + bx) k{d + 2x) 



C (l + 6x)2 (l + 6x) ' 
When k = then the mass function (1151) has the form 



(34) 



53,3/2 

"^("^ = 2C3/2(l + 6x)- 
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The expression of the mass function given in (j35p represents an energy density which is 
monotonically decreasing in the interior of the uncharged sphere and has a finite value at the 
centre x = 0. This is physically reasonable and similar mass profiles appear in the treatments 
of general relativistic equilibrium configurations of Matese and Whitman^'^, neutron star 
models of Finch and Skea^, anisotropic stellar solutions of Mak and Harko^, and the dark 
energy stars of Lobo^. The gravitational potentials for charged imperfect fluid solution 
corresponding to the line element ()34p are nonsingular at the origin. The energy density p, the 
measure of anisotropy A, and associated quantities pr and pt are also regular at the centre for 
our choice of E. Note that the quantity A vanishes at the centre and is a continuous function 
in the stellar interior. Hence the model (j27p - (|33p admits a particular case corresponding to 
the parameter values a = 1 and n = 1 which is regular and well behaved at the origin. 



4 Isotropic models 

It is possible to consider the special case of isotropic pressures with pr = pt in the uncharged 
limit for neutral matter. When k = {E = 0) equation ()32p becomes 



A 

1 r 4C(1 -a- AbCn{l + 5a)x"'i 8/? 2[4a - (n(l + a) - 4)6x'^j . ^ 

4\ x{a + bx'') (l + a)(a + 6x") ~^(l + a) {1 + a)x{a + bx"^)^ ^ 

{(3x - C(l + a)) - 6x" (C ((1 + q)(6x" - 1) +2na) + a (C(l + a)(l - 2bx^) 

+26/3x"+Ml + , ^, — - [46(n - 1)(1 + 0)2;" + (a - 1)(1 + a)(a(l + a) - 5 - q) 
x{a + bx"-) 

4(a - l)bna{l + 0)2;" 46n/?x"+i Abna[2a{n - 1) + b{na - 2)x"]x" 

(a + 6x") C {oTbx^ 

[a/3x - 6(C(1 + a) - /?x)x"]2 2 , , , 

+ (2na - (1 + a)) ((C7(l + a) - /?x)6x" - a/3x)]} . (36) 



Equation (j36p shows that the model remains anisotropic even for the uncharged case in 
general. However for particular parameter values we can show that A = in the relevant 
limit in the general solution (|27 p - (|33p . If we set 6 = and /3 = then ()36p becomes 

. C{a-l) 



4ax 



{a{l + af- [4a + (1 + a)2] } . 



From the above equation we easily observe that when a = 1 and a = 1 + (^i^y2 the measure 



of anisotropy A vanishes. When a = 1 we note from (|29p that p = since b = k = 0. 
Consequently we cannot regain an isotropic model when a = 1; to avoid vanishing energy 
densities we must have a ^ 1 when A = 0. When a = 1 + f^i^^i we obtain the expressions 



ds"^ = -Br^+<^dt^ + 



4a 

1 + 



l + a)2 



dr"^ +r'^{de'^ + s\n^ed<l)'^), (37) 



4a 

27^, Pr{=Pt) = ap, (38) 



[4a + (1 + a)2]r 
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where we have set A^D'^a'^ = B and C = 1. 

The above solution was obtained by Saslaw et al^ in their investigation of general rela- 
tivistic isothermal universes. Since p oc we may interpret (j37p as a relativistic cosmological 
metric since there is an analogy that can be made with the well known Newtonian solution 
as pointed out by Chandrasekhar— . In the Newtonian case the density p is finite in the core 
regions and decreases according to in the rest of the model. The total mass and region 
of the isothermal model are infinite. Therefore the line element (I37p may be interpreted as 
a relativistic inhomogeneous universe where the nonzero pressure balances gravity. Saslaw 
et al— point out that (j37p may be viewed as the asymptotic state of the Einstein-de Sitter 
cosmological model, in an expansion- free state as i — > oo, where the hierarchial distribution 
of matter has clustered over large scales. As t — > oo the Einstein-de Sitter model tends to 
the line element (|37p and the universe evolves into an isothermal static sphere given by the 
exact solution ([37|) with equation of state Pri= Pt) = oip. To move from the pressure- free 
Einstein-de Sitter model to the static isothermal metric requires a phase transition for conden- 
sation through clustering of galaxies. It is interesting to observe that Chaisi and Maharaj— 
and Maharaj and Chaisi— have obtained generalised anisotropic static isothermal spheres by 
utilising a known isotropic metric to produce a new anisotropic solution of the Einstein field 
equations. Govender and Govinder— have found simple nonstatic generalisations of isother- 
mal universes which describe an isothermal sphere of galaxies in quasi-hydrostatic equilibrium 
with heat dissipation driving the system to equilibrium. Hence the imperfect relativistic fluid, 
in the presence of electromagnetic fleld, with the linear equation of state pr = oip — (3, is a 
generalisation of the conventional isothermal model with a clear physical basis. 



5 Physical analysis 

We observe that the exact solution (p7|) - ([33]) may be singular at the origin in general. The 
solution should be used to describe the gravitational field of the envelope in the outer regions 
of a quark star or a dark energy star with equation of state pr = ap — (3. To avoid the 
singularity at the centre another solution is required to model the stellar core. Examples 
of core-envelope models in general relativity are provided by Thomas et al—, Tikekar and 
Thomas '^^ and Paul and Tikekar—. We observe that if a = 1 then the imperfect charged 
solution ()27p - (j33p admits the line element 

ds^ = -A^D^{1 + exp[2F{x)]dt^ + (1 + bx'')dr^ + r\d9^ + sin^ Od^'^), (39) 

where F{x) = \l + ^1 + - k{dx + x^)] . It is clear from §9^) that the 

gravitational potentials are nonsingular at the origin for all values of n. If n = 1 then there 
is no singularity in the energy density p and the model is nonsingular throughout the stellar 
interior. However if n 7^ 1 then the gravitational potential potentials (139p may continue to 
be well behaved but singularities may appear in the matter variables at the origin r = 0. 

By considering a particular example we can demonstrate graphically that the matter vari- 
ables are well behaved outside the origin. Figures [l]l5] represent the energy density, the radial 
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pressure, the tangential pressure, the electric field intensity and the measure of anisotropy, 

respectively. Note that solid lines represent uncharged matter {E = 0) and dashed lines 

include the effect of charged matter E ^ 0. To plot the graphs we choose the parameters 

n = 2,a = l,b = 40, A; = 2,d = l,a = | and /3 = 0.3569, and the stellar boundary is set 

at r = 1. From Figures [T][2] we see that both the energy density p and the radial pressure 

Pr are continuous throughout the interior, increasing from the centre to r = 0.32 and then 

decreasing. Note that the radial pressure is zero at the boundary r = 1 for the uncharged 

case E = 0. We observe from Figure [3] that the tangential press pt is continuous and well 

behaved in the interior regions. From Figured] we observe that the electric field intensity E 

is decreasing smoothly throughout the stellar interior. We can observe from Figure [5] that 

the measure of anisotropy is continuous throughout the stellar interior. The behaviour of A 

outside the centre is likely to correspond to physically realistic matter in the presence of the 

electromagnetic field. Figure 5 has a profile similar to the anisotropic boson stars studied by 

Dev and Gleiser— and the compact anisotropic relativistic spheres of Chaisi and Maharaj— . 

From the figures we can see that the effect of electric field intensity E is to produce lower 

values for p,Pr,Pt and A. 
p 
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Figure 1: Energy density. 




E = 



Figure 2: Radial pressure. 
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0.2 0.4 0.6 0.8 1.0 

Figure 4: Electric field intensity. 

We now show that the solutions generated in this paper can be used to describe realistic 
compact objects for the case n = 1 as an example from §3. In our model, when n = 1, the 
parameters b has the dimension of length"^, k has the dimension of length"^ and d has the 
dimension of length^. For simplicity, we introduce the transformations 

b = bR^, k = kR^, d = dR-"^, 

where i? is a parameter which has the dimension of a length. Under these transformations 
the energy density becomes 



+ 



2b k{d + 2y) 



l + by (1 + byY I + by 




Figure 5: Measure of anisotropy. 
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where we have set C = 1 and y 
form 



r 

7?^ 



1 



M = - 



2k 



Then the mass contained within a radius s has the 
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2 I 1 + 36 



+ k{2 - bd) 



.62 55/2 

For simphcity we set bd = 2 so that these expressions reduce to 



arctan 



P = -f;o 
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1 + by (1 + byf 



2k 
T 



M 



1 + 6^ 



2k 
36 



s 

R^' 



(40) 
(41) 



which are simple forms. It is now easy to calculate the density and mass for particular 
parameter values from (|40p and (j42p . For example, when s = 7.07km, R = 1km, 6 = 0.03 and 
k = 0.00045, we obtain the mass Me=o = 1.436M0 for uncharged matter MejLq = O.24OM0 
for charged matter. Note that the value of the mass for uncharged matter is consistent with 
the strange star models previously found by Sharma and Maharaj— and Dey et al^^. We 
have shown that the inclusion of the electromagnetic field affects the value for the mass M. 
It is also possible to to relate our results to other treatments. If we set s = 9.46km, R = 1km, 
6 = 0.35 and k = 0.00045, then we obtain the mass Me=o = 3.1O3M0 for uncharged matter 
and Me^o = 2.858M0 for charged matter. These values for the mass are similar to charged 
quark stellar models generated by Mak and Harko^ which describe a unique static charged 
configuration of quark matter admitting a one-parameter group of conformal symmetries. 
We have shown that the presence of anisotropy and charge in the matter distribution yields 
masses which are consistent with other investigations. 



6 Discussion 



We have investigated the Einstein-Maxwell system of field equations with a strange mat- 
ter equation of state for anisotropic matter distributions when the electromagnetic field is 
present. A new class of exact solutions was found to this system of nonlinear equations. A 
particular case is nonsingular at the centre of the star and the gravitational potentials and 
matter variables are well behaved. When the anisotropy vanishes we regain the general rela- 
tivistic isothermal universe of Saslaw— . This is an extension of the conventional Newtonian 
isothermal universe with density p oc r^'^. The plots of the matter variables indicate that the 
solution may be used to model quark stars, at least in the envelop if singularities are present 
at the origin. We calculated the mass in a special case (a = l,n = 1) and showed that 
this value is consistent with strange matter distributions of Dey et al-^, Mak and Harko^ 
and Sharma and Maharaj—. Ours is a particular solution of the Einstein-Maxwell system 
with the nice feature of containing the isothermal universe in the isotropic limit. A more 
comprehensive study of other possible solutions admitted, with the strange matter equation 
of state, is likely to produce other interesting results. 
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We now comment on some physical aspects of the solutions found which are of interest. 
Firstly, there are a number of free parameters in the solution which may be determined by 
imposing boundary conditions. The solutions may be connected to the Reissner-Nordstrom 
exterior spacetime 

rf^' = - (l - ^ + J) + (l - ^ + 5) ' + (dO' + sin' ^^'^') (42) 

with the interior spacetime (j39p across the boundary r = s where M and Q represent the 
total mass and the charge of the star, respectively. This gives the conditions 

2 /If r)2 

1- — + ^ = ^'^M[l + KC*^)"]"exp[2F(C7s2)]| 

Therefore the continuity of the metric coefficients across the boundary r = s is maintained. 
The number of free parameters {A, C, D, b and n) ensures that these necessary conditions are 
satisfied. Secondly, from the graphs generated for the energy density and the pressure we 
observe that these quantities are negative in regions close to the centre. We believe that this 
feature arises because of the strange equation of state that we have imposed on the model. 
To avoid negative matter variables we would need to use another solution to describe the core 
with our solution serving as the envelope. This is the situation that arises in the gravastar 
model of Mazur and Mottola where the singular core is replaced by a de Sitter condensate 
through a phase transition near the location of the event horizon. This idea has been utilised 
by Lobo^^ to demonstrate that dark energy stars are stable. Thirdly, the example of this 
paper and other exact solutions found for the Einstein-Maxwell system suggests that the 
negativity of the energy density and the pressure close to the centre may be a generic feature 
of charged stellar models. It appears that the existence of the electric field always makes the 
energy density negative close to the centre of the sphere. This suggests the conjecture that 
a charged fluid cannot satisfy the energy conditions near the stellar centre is possible. This 
will be the topic of a separate investigation. 
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